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ABSTRACT 

We present the two-dimensional (2D) ionization structure of self-similar mag- 
netohydrodynamic (MHD) winds off accretion disks around irradiated by a cen- 
tral X-ray point source. Based on earlier observational clues and theoretical 
arguments, we focus our attention on a subset of these winds, namely those with 
radial density dependence n(r) oc 1/r (r is the spherical radial coordinate). We 
employ the photoionization code XSTAR to compute the ionic abundances of a 
large number of ions of different elements and then compile their line-of-sight 
(LOS) absorption columns. We focus our attention on the distribution of the 
column density of the various ions as a function of the ionization parameter £ (or 
equivalently r) and the angle 9. Particular attention is paid to the absorption 
measure distribution (AMD), namely their Hydrogen-equivalent column per log- 
arithmic £ interval, dNff/dlogl;, which provides a measure of the winds' radial 
density profiles. For the chosen density profile n(r) oc 1/r the AMD is found to 
be independent of £, in good agreement with its behavior inferred from the X-ray 
spectra of several active galactic nuclei (AGNs). For the specific wind structure 
and X-ray spectrum we also compute detailed absorption line profiles for a num- 
ber of ions to obtain their LOS velocities, v ~ 100 — 300 km s _1 (at log£ ~ 2 — 3) 
for Fexvil and v ~ 1,000 - 4,000 km s^ 1 (at log£ ~ 4-5) for Fexxv, in 
good agreement with the observation. Our models describe the X-ray absorption 
properties of these winds with only two parameters, namely the mass-accretion 
rate m and LOS angle 9. The probability of obscuration of the X-ray ionizing 
source in these winds decreases with increasing m and increases steeply with the 
LOS inclination angle 9. As such, we concur with previous authors that these 
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wind configurations, viewed globally, incorporate all the requisite properties of 
the parsec scale "torii" invoked in AGN unification schemes. We indicate that a 
combination of the AMD and absorption line profile observations can uniquely 
determine these model parameters and their bearing on AGN population demo- 
graphics. 

Subject headings: accretion, accretion disks — galaxies: active — methods: nu- 
merical — quasars: absorption lines — X-rays: galaxies 



1. Introduction 

The issue of accretion as the power sour ce behind Active Galactic Nuclei (AGNs) was 



decided as early as 1969 (jLynden-Belll Il969l ). Since then, a 40- year effort to unravel the 



physics underlying this process has produced a great body of observational and theoretical 
work covering a host of the AGN properties and phenomena. However, due to the small 
angular size of a black hole horizon, whether at the galactic or extragalactic setting, these 
studies were effectively performed mainly in the spectral (but also the time) domain; the 
spatial structure of accretion flows, but those of scale of many parsecs, was then inevitably 
delegated to models, whose validity was determined by their ability to reproduce the spectral 
and timing observations with sufficient accuracy. At the risk of oversimplifyin g the issue, the 



spectr al properties of AGNs were basically outlined (among other works) in ISanders et al. 



( 119891 ) who showed that the Spectral Energy Distribution (SED) of AGNs includes three 
broad components in the IR, optical-UV and X-ray bands, with roughly equal energy per 
decade, with the radio comprising only a small fraction of their bolometric luminosity, of 
order ~ 10" 6 in the Radio-Quiet (RQ) AGNs and ~ 1(T 3 in Radio-Loud (RL) AGNs. 

Considering that accretion produces most of its power in the last 10 or so Schwarzschild 
radii, r$, outside the black hole horizon, barring an inherently non-thermal emission process 
(as it turns out to be the case with blazars), it is strange that the observed SEDs exhibit 
roughly constant luminosity per decade. Interestingly, the simplest of assumptions, namely 
radiation of the AGN bolometric luminosity in black body form by an object (most likely 
in the form of a cold disk of T ~ 10 4 — 10 6 K) of size a few of a supermassive black hole 
(M ~ 10 6 — 10 s M & where M & is the solar mass), implies peak emission at UV frequencies, 
consistent with the observationally identified UV feature referred to above, known as the 
Big Blue Bump (BBB); as such, this feature was modeled with emission by a geometrically- 
thin, optically-thick accretion disk extending usually to the Innermost Stable Circular Orbit 
(ISCO) of the black hole. Of the other two distinct components the X-rays are generally 
attributed to emission by an optically-thin, hot (T ~ 10 8 — 10 9 K) corona overlying this disk 
and heated by the action of magn etic fields that thread the geometrically-thin disk, in close 



analogy to the Solar Corona (e.g. lHaardt fc Maraschilll99ll : iKawanaka. Kato. fc Mineshige 
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20081 ). Finally, the IR emission was accounted for as the result of reprocessing the O- 
UV radiation by a geometrically thick (scale height h roughly equal to the local radius r, 
i.e. h ~ r), cool, molecular torus at large distances (> 1 pc) from the central engine [an 
arrangement that presents a bit of a problem, since the temperature of this torus (~ 10 — 100 
K) is much less than the virial temperature of the gas at that distance and would lead to a 
configuration with h <C r; see Krolik & Begelman (1988). The region between the torus and 
the accretion disk was thought to be occupied by clouds with velocity widths ranging from 
tens of thousands to several hundreds of kilometers per second, producing the observed line 
ra diation, with the e nsemb le of these components encapsulated in the well known cartoon 
of lUrry &: Padovanil (119951 ). 



While in most models these components are generally thought of as i ndependent, Prin- 
cipal Com ponent Analysi s (PCA) of multiwavelength AGN features by iBoroson fc Green 
fll992h and lBorosonl f|2Q02h indicated that most of the variance in the measured optical emis- 
sion line properties and a broad range of continuum features (radio, optical, X-ray) of AGNs 
was contained in two sets of correlations, eigenvectors of the correlation matrix. This analy- 
sis suggested that the properties of these features are not independent (despite their diverse 
locations and emission processes) but are well correlated by two single underlying physical 
parameters through the physics of accretion and the conversion of its power to radiation and 
outflows. On the other hand, general theoretical arguments suggesting, e.g. a tight correla- 
tion between variations in the X-ray and O-UV compon ents, given their implied proximity, 
were not confirmed by observations (iNandra et al.lll998l ). 



The advent of X-ray instrumentation and in particular X-ray spectroscopy, established 
the presence of spectral components apparently present across the entire range of compact 
object masses from galactic black hole candidates (GBHCs) to luminous AGNs, namely out- 
flowing X-ray absorbing matter in our line of sight (LOS) manifested as complex (blueshifted) 
absorption features of hydrogen equivalent column density Njj ~ 10 21 — 10 23 cm -2 . The pres- 
ence of this absorbing material, often referr ed to as Warm Absorber, w as first established 
in ASCA observa tions of many AGNs (e.g. iReynolds and Fabianl Il995l ) and GBHCs (e.g. 
Miller et al.l 120061 ) which provided clear evidence of highly ionized oxygen O VII (0.739 keV) 
and/or Ovill (0.871 keV) along the LOS. More recently, grating spectra of high-resolving 
power obtained by Chandra and XMM-Newton have enabled the study of these absorption 
features in greater detail leading to the conclusion that they are pre sent in a large fraction of 
AGN s and span a wide range (~ 10 5 ) in ionization parameter (e.g. iHolczer. Behar fc Kaspi 
20071 . hereafter HB K07). Some sources appear to contain low and/or high velocity out- 



flows (e.g. HBK07; IHolczer. Behar fc 



relativistic outflow (e.g. IChartas et al. 



Arav 



2009 



2009) while a handful of objects exhibit a tr ans- 



Pounds fc Reeved BOOalReeves et al.ll2009h . In 



addition in both AGNs and GBHCs the observed outflows from at least a number of sources 
have been suggested to be accelerated magnetically (rather than by another process) off an 
accretion disk (e.g., see Miller et al. 2006, 2008 for GRO J1655-40; Kraemer et al. 2005 and 
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Cr enshaw fc Kraemer 2007 for NGC 4151). A brief review of luminous AGN winds is found 



in 



Brandt et al.l ( 120091 ). These facts imply the presence of gas covering a large fraction of 



the solid angle and distributed over a large range of radii, perhaps the entire range between 
the X-ray source and the torus, indicating AGNs to be multiscale and multiwavelength ob- 
jects rather than a class with properties determined only by the power released by accretion 
in the black hole vicinity. This last point has gained further support with the discovery 
of the correlation between the accreting black hole mass and the velocity dispersion of the 
surrounding stellar population ( IFerrarese fc Merritll2000l ). 



From the theoretical point of view, most early AGN treatments were focused on the 
structure of the innermost regions of the accretion flow and as such they limited themselves 
to the study of the conditions in a rather narrow range of radii. Nonetheless, several treat- 
ments of magnetohydrodynamic (MHD) accretion disk winds, aiming to account for the 
observed AGN outflows, employed self-similar solutions, which naturally span a large range 
in radius (e.g. Blandford & Payne 1982, hereafter BP8 2; Contopoulos & ; Love lace 1994, 
hereafter CL94). Using the structure of these solutions, iKongl &: Kartjd (119941 . hereafter 
KK94) proposed that the so-called molecular torus of the AGN unification scheme is a dy- 
namical rather than a static object, namely an MHD wind of the type suggested in the above 
works. KK94 also contented that agreement between model and observation demanded a 
rather specific type of wind, a particular case of those discussed in CL94, one that we also 
concentrate on in the present work. At the same time, they noted that the 2D geometry of 
these winds (with low column at inclination angles 9 ~ 0°, i.e. along the wind axis, and high 
column at 9 ~ 90°) provided, in addition, a natural framework for the unification scheme of 
type I and II Seyfert galaxies proposed by Antonucci & Miller (1985), thereby linking the 
large scale AGN spectral classification to the dynamics of AGN accretion/outflows. 



S elf-simi lar solutions of the accretion flow equations were also provided by lNarayan fc Yi 



( 11994J . Il995al ). who considered the structure of hot accretion flows in the regime of low accre- 
tion rate, i.e. for accretion rates less than the Eddington value, m = M/Me ^ 1. Because 
for these low rates the accretion time scale r acc is shorter than the gas cooling time r coo i, in 
fact r acc ~ mr coo i, only a fraction fa of the energy released in the dissipation of the plasma 
kinetic energy is radiated away with the remainder advected into the black hole (hence the 
term Advection Dominated Accretion Flows or ADAFs); as a result, the accretion luminos- 
ity is then proportional to fn? (rather than fn). The high temperature of the inner ADAFs, 
provides naturally for the hot component required to produce the observed X-ray emis- 
sion, which, rather than being an in dependent component, within ADAFs is incorporated in 
the dyna mics of accretion itself ( see iNaravan fc McClintockl 120081 . for recent review). More 
recently, Blandford Sz Begefmanl ( 19991 ). elaborated further on ADAFs elucidating, among 
others, the re ason for which the Bernoull i integral of ADAFs is positive (a point that had 
been noted by lNarayan fc Yilll994l . Il995al ). that being that the viscous stresses transfer out- 
ward, in addition to angular momentum, also mechanical energy. They then proposed that 
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this excess energy can be carried away in the form of wind over a large range of radii, thereby 
leading to configurations similar to those obtained in the more detailed studies of BP82 and 
CL94, while at the same time maintaining some of the general properties of ADAFs. The 
combination of radiatively inefficient flows with the simultaneous presence of winds led to 
the nomenclature Advection Dominated Inflow-Outflow Solutions (ADIOS). Following these 
pioneering works, many attempts have been made in recent years to reproduce and explain 
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Motivated by these considerations we are taking a closer look at the issue of Warm 
Absorbers (X-ray absorbing medium) in terms of specific disk-wind models, namely those of 
CL94. As it will be discussed in more detail later on, these solutions are characterized by 
a parameter q, which determines the distribution of axial current in the wind as a function 
of the radius; this parameter, in conjunction with the MHD conservation laws and the 
radial force balance, determines both the radial dependence of the matter density in the 
wind and the radial dep endence of t h e toro idal magnetic field component B^. Taking our 
lead from the works of iBehar et al.l ( 120031 ) and HBK07, in the present work we restrict 
our attention to wind models with q = 1, i.e. the value used also by KK94. This is in 
a sense a critical value of this parameter because, as discussed in CL94, it leads to winds 
with toroidal magnetic field B$ oc 1/r and therefore magnetic field energy on the disk that 
diverges only logarithmically at small and large r and as such it can be considered as a 
"minimum magnetic energy" configuration (KK94); values of q ^ 1 lead to configurations 
with power-law divergence either at small or large R. In addition, the choice q = 1 leads to a 
wind with radial density profile n(r) oc r -1 , i.e. a profile with equal column per logarithmic 
radial interval (and normalization that depends on the polar angle 9 and the mass-accretion 
rate m), a fact that is in general agreement with the profiles implied by the ionized absorbers 
da ta analysis to date (HBK07) . It is wort h noting that a similar densi t y profi le was invoked 
by iKazanas. Hua &: Titarchukl (119961 ) and iHua. Kazanas &: Titarchukl (119971 1 to explain the 
Fourier frequency dependence of the soft-hard X-ray lags observed in GBHCs and in AGNs 
(Papadakis et al. 2002). 

With the wind 2D (r and 6) density field as a function of radius and angle provided by 
the models of CL94 (i.e. ignoring at present effects that can be attributed to other agents 
that could affect the wind structure) and the assumption of a point X-ray source of a given 
spectrum at the origin, one can compute their 2D ionization structure and in particular the 
local column density of specific ionic species as a function of radius and angle. One can also 
compute the integrated column of each such ion, quantities that can be directly compared 
to observation (HBK07). This we do in the present work. While a number of issues remain 
open or are sidestepped in the present self-similar wind models, this is the price to pay in 
order to limit the number of free parameters to just two (the wind mass flux and the observer 
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inclination angle), a fact that allows the study of the properties of warm absorbers within 
the global perspective of AGN unification. 

Our paper is structured as follows: in §2 we review the MHD wind equations and struc- 
ture as given in CL94 and make a connection of these winds with flows of the ADAF/ADIOS 
type by relating the wind density normalization to the ionizing luminosity produced by the 
idealized point source at the coordinate origin. In §3 we provide the general 2D ionization 
structure, the local column of specific ions as a function of radius as well as the total ionic 
column of a given charge state and for each element, as well as the absorption line profiles 
of selected transitions. In §4 the results are compared to observation, the properties and 
limitation of the present treatment are discussed and a course for future work is charted. 



2. Description of the Model 



Our wind model consists primarily of two parts: (i) The outflow/wind structure origi- 
nating from a geometrically-thin accretion disk as described by CL94 and (ii) The compu- 
tation of the ionization of the wind plasma under local heating-cooling and ionizat ion equi- 
librium, by employing the photoinization code XSTAR (IKallman fc Bautistal 120011 ). These 
two steps are self-consistentl y coupled in a scheme similar to some previous work in a 
slightly different context (e.g . ISchurch &: Pond 120071 : iDorodnitsyn. Kallman fc Prog 



Schurch. Done fc Progall2009i ). All calculations are performed under steady-state, axisym- 
metric conditions as shall be described in detail in the following sections. 



2.1. Self-similar MHD Wind from Accretion Disks 

While the detailed formalism of self-similar MHD wind solutions is found in CL94, let 
us here briefly summarize some of their characteristics. The basic equations of this problem 
are those of steady-state, nonrelativistic, ideal MHD that include gravity and gas pressure, 
namely 



V ■ (pv) = 


(mass conservation) , 


(1) 


„ 4:71 

V x B = — J 
c 


(Ampere's law) , 


(2) 


E+ - x B = 
c 


(ideal MHD) , 


(3) 


V x E = 


(Faraday's law) , 


(4) 


p(v ■ V)v = -Vp - pV$ g + i(J x B) 


(momentum conservation) , 


(5) 


VB = 




(6) 
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Here, the plasma wind velocity v is assumed to be frozen into the global magnetic field B; 
E is the electric field; J is the electric current density; and c is the speed of light. The gas 
pressure p and mass density p of the wind are related through an equation of state 

p = Kp r , (7) 

where K and Y are the adiabatic constant and adiabatic index respectively. The assumption 
of axial symmetry allows one to separate poloidal quantities from toroidal ones as B = 
Bp + B^ and v = v p + v^, with 

Bp = -^-Wx0, (8) 
r sin 9 

where, \&(r, 9) is the magnetic flux function. We will henceforth work in spherical coordi- 
nates (r, 9, 4>) [note that CL94 used cylindrical ones (R, 0, Z)\. Magnetic field lines lie along 
surfaces of constant value of Furthermore, under conditions of steady-state, the following 
quantities are conserved along field lines characterized by the value of ^ 

4^ ee , (9) 

— ?-i (v+ - -^-bA ee n(tt) , (io) 

rsmv \ Anp J 
rsm9(B (j> -v^F)=H(^) , (11) 

/ (dp/p) + \v 2 + $ g - r sin 9v 4> Q = J(tt) , (12) 

J <&=const ^ 

(r — 1) 1 \n(K) ee S(ty) , (13) 

where F(^) is the ratio of mass to magnetic fluxes, is the angular velocity of mag- 

netic field lines, H(fy) corresponds to the specific angular momentum (it includes magnetic 
torques), J(^) is the specific energy (Bernoulli) integral, and S(^f) is the specific entropy. 
Also, Jcb is Boltzmann's constant, and $ 9 is the Newtonian gravitational potential. The 
wind velocity field v is related to the magnetic field B as 

v = rsin0Q0+— B. (14) 
4irp 

CL94 first showed that one can obtain general self-similar solutions of the above system of 
equations of the form 

B(r,9) ee {r/r f- 2 h{9)B , (15) 

v(r,0) ee (r/r o yV\(9)v , (16) 

p(r,0) ee (r/r ) 2 ^(^)5 o 2 , (17) 

pM) ee (r/r ) 2 '^(^ \- 2 . (18) 
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Here, r Q is the launch radius of a characteristic field/flow line at the base of the wind at 
9 = 90° while B and v are respectively the magnitude of the vertical components of the 
magnetic field and the initial rotational wind velocity at the same foot point; b(9) and 
v(9) are dimensionless with 6^(90°) = v^(90°) = 1 while V(9) and 11(9) denote the angular 
dependence of the pressure and wind density, respectively. The exponent q is a free parameter 
that controls the radial scaling of the magnetic field. One sees directly that q discriminates 
current-carrying wind configurations (when q > 1, rB^ grows with r; the wind carries a 
net axial current) from zero-total-current ones (when q < 1, rB^ — y for r — > oo; the 
wind carries a singular axial current and its corresponding return current; see CL94 and 
references) . Obviously, 

*(r,0) = (r/r o )'V(0)* o , (19) 

where ip(9) is dimensionless, ip(90°) = 1, and = B o rlb z (90°)q~ 1 is the poloidal magnetic 
flux through radius r Q . The above scalings allow us to obtain the geometry of the poloidal 
field/flow lines (characterized by the value of ty) as 

r|*(0) = r (^/^ ) l ^(9Y 1 l q . (20) 

Note that the function ip(9) characterizes the shape of all field/flow lines. That is, all 
field/flow lines are self-similar to the 'characteristic' field/flow line that originates at r = r Q 
and 9 = 90° on the disk. Self-similarity further implies that 



F(¥) = 




'-'^FoBoV- 1 , 


(21) 


n(*) = 


(*/*o) 




(22) 


Hm = 




l - lt(l H B r , 


(23) 


m = 






(24) 


1) _1 [{(2- 


- 4/g) - 


r(2-3/g)}ln(^ ) + ln/C] , 


(25) 



where K = V(9)/1Z(9) r is a dimensionless adiabatic constant. The adiabatic speed of sound 
scales as 



dp 



dp 



1/2 



c s (r,9) ee ( ^ ) = (*/¥ )-V2* /Cr^(0) i " i ' t; G = (tf/^r^c^K . (26) 



^=const. 



r-i 



1/2 



Finally, one can combine the dimensionless forms of equations (I9l-(l25l) that involve the 
dimensionless quantities (v, b, V, 1Z, c s ) as functions of 9, and thus express the poloidal pro- 
jection of the momentum balance equation (jSJ) [the so called generalized Grad-Shafranov 
(GS) equation] in the form of a second-order ordinary differential equation for ip(9) as 



f(9) = f(9^^') 



(27) 
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This equation (too complicated to write down explicitly in this paper) is integrated numer- 
ically from the surface of the disk (taken for simplicity to lie at 9 = 90°) to axial infinity 
(9 = 0°). The initial conditions are ^(90°) = 1 and ip'(90°), essentially a free-parameter 
that depends on the details of the internal disk structure, which, for a given value of H Q , is 
adjusted in order for the wind to cross the Alfven point (the point at which the the wind 
poloidal speed v p becomes equal to the local Alfven speed - see CL94 for a detailed descrip- 
tion) continuously and smoothly. According to equation (1201 . the form of the solution tp(9) 
characterizes the poloidal shape of the field lines, and in particular, the initial condition 
?/>'(90°) characterizes the angle the poloidal field/flow lines make with the disk at their foot 
points [readers can easily verify that this angle is equal to tan~ 1 {g/-?//(90 )}]. 

To numerically solve the above wind equation (l27jl we assume that the plasma is initially 
in Keplerian rotation v Q = 1^(90°) = v k at the disk surface (9 = 90°) with only a small out-of- 
the-disk velocity component (v r ,Vg <C at 90°). As a fiducial solution we fix the following 
set of parameters as described in CL94; q — 1, F a — 0.15, H a = —1.584 with kinematic 
variables ^(90°) = 1, v e (90°) = 0.01, /C = 0.01 and T = 5/3. With these initial conditions, 
the remaining conserved quantities are calculated as fl D = 1 — vg(90°)(F o + H a ) = 0.983, and 
J = v 2 9 (90°) [1 + (^'(90°)) 2 ]/2 - 1/2 - n o ~ -1.483 (see CL94 for details). The only freedom 
left is the value of ip'(90°) that is obtained iteratively so that the wind passes smoothly 
through the Alfven surface. This important physical constraint yields ^'(90°) = 0.98091. 
Note that for sufficiently small F a (a magnetically dominated flow) its velocity is nearly 
Keplerian, setting the value of fl a , while passing through the Alfven point fixes H a and 
effectively the remaining constants of the flow, with q as the one meaningful free parameter. 
The obtained wind kinematics are shown in Figure [Tk^ and the magnetic field components 
in Figure [lb. As discussed in CL94, different values of the exponent q allow for distinct 
streamline geometries in the poloidal plane, which in turn is associated with other observable 
quantities of the wind (such as the density profile; see following sections). As is initially set 
up, this wind is launched from a thin Keplerian disk with only a small out-of-the-disk surface 
velocity component. The magnetic field there is dragged by the rotating plasma (v^ > 0) 
resulting thus in a negative toroidal field component (B^ < 0). As the wind leaves the disk, 
magnetic torques act on the outflowing plasma and magnetic energy is efficiently converted 
into kinetic energy of the wind along the symmetry axis. The magnetic field also plays an 
important role in collimating the wind at large distances. In this particular case the wind is 
found to end up with terminal velocity roughly a few times the initial rotational speed and 
remains below fast magnetosonic speed. 

The astute reader may note that X-ray heating is likely to produce thermal velocities 
inconsistent with the self-similarity of the sound speed invoked above (i.e. c 2 s oc 1/x; see 
Fig. HJ). We will discuss this in §3.1. 

Knowing the wind streamline described by equation ( 1271) one can separately write down 
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Fig. 1. — MHD wind solution of the GS-equation ff27|) in the poloidal plane: (a) poloidal v p 
(solid) and toroidal (dashed) velocities (in units of v = vk) as a function of axial distance 
Z = (r/r D ) cos 9 and (b) poloidal B p (solid) and toroidal \B^\ (dashed) magnetic fields (in 
units of B ) also as a function of Z along the characteristic field/flow line originating at 
r = r on the disk for the parameters described in the text. Note that we show the absolute 
value of B^(< 0). Vertical dotted lines denote the axial location of the Alfven point. 



radial and angular dependence of the wind number density n(r, 9) as 

n(r, 9) = = n x 2q - 3 Af(9) , (28) 

where n = B 2 TZ(90°) / '(/j,m p v^) is the number density normalization at the initial character- 
istic launching radius r Q , N{9) = TZ(9)/7l(90 o ) is the normalized angular dependence of the 
wind number density, and x = r/r Q is the non-dimensional radial coordinate; m p is the pro- 
ton mass and \i = 1.26 the mean molecular weight of the wind. In what follows, we will take 
the characteristic radius to be equal to the Schwarzschild radius, namely r Q = r$, taken to be 
the inner edge of the accretion flow with r$ — 3 x 10 5 M cm where M = M/M & is the mass 
of the black hole in units of solar mass M & . It is well known that free-fall spherical accretion 
at the Eddington rate (m = 1) produces Thomson opacity one (tt — 1) at the Schwarzschild 
radius (r ~ rg). Using this scaling we can express our wind density normalization at the 
inner edge of the disk as 

n Q = , 29 

2a T r s 

where, i]w is the ratio of of the mass-outflow rate in the wind to the mass-accretion rate 
m, assumed here to be r]y/ ~ 1, and gt is the Thomson cross-section. It is important to 
note here that because the mass flux in these winds depends in general on the radius, a 
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Fig. 2. — Wind properties along a streamline (of \I/ = \l/ ) as a function of LOS angle 9: 
(a) density profile n{0)^=^ o (normalized to unity) given by equation ( 1271) . (b) LOS wind 
velocity t>i os along the same characteristic field/flow line (in units of v Q = v K at the base of 
the wind). A dot denotes the position of the Alfven point. The parameters are the same as 
in Figure [TJ 

normalized parameter used throughout this work, m, always refers to the mass flux at the 
innermost value of the flow radius, i.e. at x ~ 1. 

Figure exhibits a LOS angle (9) dependence of the wind density profile n(9)^ = ^ o 
from the pole (9 = 0°) to the equator (9 = 90°) along a characteristic stream line of \1/ = \l/ , 
normalized to unity at its highest value at 9 = 90°, as obtained from the self-similar solution 
of Figure [TJ This distribution illustrates the significant change of the LOS density (and also 
column density) with inclination angle (over three decades between 5° < 9 < 90°). The 
LOS wind velocity t>i os (in units of v Q = vk) along a characteristic streamline is shown in 
Figure [2b. We see that v \ os increases as the LOS angle 9 becomes smaller (toward the polar 
region) because the wind is magnetocentrifugally accelerated along the way, and furthermore, 
the direction of t>i os becomes more and more parallel to the total wind velocity v as well [e.g. 



In terms of the scalings of our wind model, the total equivalent hydrogen column density 
Nh of wind over the LOS length scale of Ar is given by 



where the wind is considered to extend from an inner radius r = r- m = r Q ~ rs over a length 
scale of Ar = Ax-r$ along the LOS. One should note that the q = 1 configurations have the 
interesting property of equal column density per decade of radius, i.e. dNn/dlogr = const oc 



v los (30°) ~ 2*; los (60 )]. 




(30) 
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J\f(9) [as can be seen from Fig. [2^, n(9)^ = ^ ~ e -{Q-w )/i3 j s an exce llent fit to the angular 
dependence of the density for 9 <L 10° where 9 is measured in degrees]. For rh ~ 1 and 
9 ~ 90° these winds have a column density of ~ 10 24 cm -2 per decade of radius. As a result, 
configurations that extend over a number of decades in radius can be quite optically-thick 
in these directions even for rh <C 1. 



2.2. The Wind Ionization Structure 

With the poloidal field configuration obtained by solving equation fl27|) and the cor- 
responding wind kinematics, we now are in a position to consider its ionization structure 
assuming the presence of a point-like X-ray source at the origin. To this end we employ the 
photoionization code XSTAR version 2.1kn9, which solves simultaneously for the ionization 
balance and thermal equilibrium of the locally illuminated wind. 

The implementation of the above code involves a number of issues and leads to certain 
wind observables, amongst which are the following: (i) The ionization is locally determined 
by the photon flux per electron flux; a proxy for that is the ionization parameter £ = Lj (nr 2 ), 
with L being the luminosity of the ionizing source; (ii) The computation of the ionizing flux 
requires the transfer of radiation from the source through the LOS wind, given that part of 
it can be removed by scattering and absorption; (in) Given that our wind models produce, 
in addition to the density n(r,9), also the wind velocity field v(r, 9) we are then able to 
compute also detailed absorption line profiles; (iv) The broad range of the wind density with 
radius leads to ions with very diverse ionizatio n states, who se distribution per logarithmic 



£ interval, the AMD, discussed in HBK07 and iBeharl (120091 ) . can be used to determine the 



wind radial profile. We will discuss each one of these issues in detail below. 



2.2.1. The Ionization Parameter 

With the density profile at hand, the wind ionization structure, i.e. the local value of 
£, requires also an expression for the X-ray luminosity. Since the models we present are 
global, i.e. they cover a large number of decades in radius, it makes more sense that we 
do not detach the ionizing luminosity from the global wind dynamics. In order to do so we 
need a prescription that connects the wind mass flux rate rh w and the ionizing luminosity 
L. The simplest assumption is that Locm = rh^ = rh w , i.e. to assume that the wind mass 
rate is proportional to the accretion rate onto the black hole (we assume they are equal), 
resulting in luminosity L proportional rh. If this were indeed true, then the ionization 
structure of the wind would be independent of the accretion rate, since rh would drop out 
of the expression for the ionization parameter £ = L/nr 2 [c.f. eqn. ( 129]) ]. However, this 
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is contrary to observations (e.g. lUeda et al.l 120031 ; iTueller et al.ll2008l ). which indicate that 



the probability of source obscuration decreases with increasing source luminos ity (one should 



note, however, that obscuration may be also due to gas in galactic mergers; e.g. lHopkins et al. 



20051 ). To incorporate this fact, we provisionally adopt the prescription L oc m 2 , even though 
any power near 2 or greater would qualitatively suffice. We note that L oc fn 2 p rescription 



would be appropriate for flows such as those discussed in iNarayan fc Yi fll995bh (provided 



that rh < a 2 ; a is the usual disk viscosity parameter) with L oc m otherwise; considering 
that the latter case is valid for only a small range in m (0.1 < m < 1) than the former, we 
adopt the former, bearing in mind its limits of applicability. Expressing the luminosity in 
terms of its Eddington value we can write 

L ~ em 2 L M = 2vr em 2 HJ^l. , (31) 

where Lq ~ 1.28 x 10 38 erg s _1 is the Eddington luminosity of a one solar mass an accreting 
object and e is the efficiency of conversion of mass into radiation for m — 1. Given the 
ionizing luminosity L, and assuming the source to be point-like, located at the coordinate 
origin, the local value of the ionization parameter (ratio of ionizing flux per electron flux) is 

n(r, U)r z 

Using the luminosity, density and radius scalings given above this reads as 

t(r 6) 4?re m P c3rh ~ e 3xl0 8 m 

qy ' ; N{6)r) W x 2 *- 1 Af(9)r] W x 2 ^ 1 ' 1 ; 

Equation fl33|) thus provides the explicit dependence of the ionization parameter of the cor- 
responding wind as a function of the accretion rate m, the (normalized) distance from the 
source x and the LOS angle 9 for steady-state accretion at a rate fn into the black hole 
and actual mass outflow rate (at radius x ~ 1) rjwiTi in the wind. Figure [3] shows poloidal 
color maps of (a) the wind density n(r,8) and (b) the ionization parameter £(r, #) for the 
parameters used to produce Figure [1] along with a number of contour curves (dotted curves 
with values) and the flow streamlines (solid curves). The diagonal dashed straight lines 
denote the self-similar surface of the Alfven point. Here, we have assumed Mbh — 10 6 M Q 
and L = 3.3 x 10 42 erg s -1 , which for e ~ 0.2 implies fn ~ 0.1. As seen, we obtain a 
non-recollimating, non-oscillating wind which asymptotically behaves as r oc (cot O) 1 ^ sec 9 
where ( = const. ~ 1/2 as r — > oo [or in cylindrical coordinates (Z,R) R oc Z 1 ^^; see CL94]. 
We would like to stress that the field line geometry depends on the wind's conserved quan- 
tities in the model. As a consequence the obtained wind properties (shown here in Figs. [1] 
to [3]) can change quite a bit. 

In the context of our accretion-disk wind scenario the outflows are necessarily more 
opaque near the launching disk surface {6 = 30°) at smaller radius as expressed in equa- 
tion ([28]) . Density profile along a magnetic field follows exactly the one shown in Figure^. 
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Fig. 3. — Poloidal distribution of (a) the wind density (in units of cm -3 ) log n(r, 9) and 
(b) the ionization parameter (in units of erg cm s" 1 ) log£(r, 9) (dotted curves; the solid 
lines denote the poloidal magnetic field/velocity stream lines; the straight dashed line is the 
Alfven surface) in a 10 8 r o x 10 8 r o region for an ionizing luminosity L = 3.3 x 10 42 , erg s _1 ; 
this calculation of £ ignores the reduction of L at large r due to absorption. The values of 
the corresponding quantities are noted on each dashed curve. The parameters are the same 
as those of Figure [TJ [See electronic edition of the Journal for a color version of this figure.} 



There are several points to note about equation (|33|) and Figure [3b: (i) If distances are 
scaled by the Schwarzschild radius, rg, and the accretion rate by its Eddington value, M^dd, 
the ionization parameter, £(r, 9), depends primarily on the normalized accretion rate m and 
H{9)\ as a result, our analysis and the corresponding ionization structure of the wind, are 
applicable to accreting black holes of any size. Note that, had one assumed Locm, the ex- 
pression for £ would be independent even of m; we believe that this would be in disagreement 
with observation; (ii) The ionization parameter £ is much higher near the central engine (e.g. 
within a few hundreds to thousands Schwarzschild radii) with the wind highly ionized in this 
region and for a given r it drops precipitously with 9; {Hi) For a given 9, the value of £ drops 
with distance (with £ oc 1/r for our favorite value of q ~ 1); (iv) £ increases linearly with 
increasing accretion rate m, implying decreasing sourc e obscuration for h igh values of rh as 
seen in equation fl30l) . in agreement with observation (ITueller et al.ll2008l ); (v) The wind is 
fully ionized in regions sufficiently close to the central engine thus no contribution to the 
absorption lines in the spectrum. 
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One should note that, as written, the expression of equation (133|) is a local one, i.e. it 
does not take into account neither the absorption of the X-ray radiation in its propagation 
along the LOS wind nor the corresponding re-emission, both of which would modify the 
local ionizing spectrum, and thus also the global ionization structure shown in Figure [3b 
here. These effects are small at small values of the radius (logx < 3 and for relatively hard 
spectra a ~ 1 — 2), where the ionization of the gas is high and its absorption depth small 
at all energies; they become more significant at larger radii and they are taken into detailed 
consideration in our model calculations as discussed below. 

It should also be cautioned at this point, that while £,{r,8) is the primary parameter 
that determines the ionization state of the X-ray illuminated gas, the latter depends also on 
the spectrum of the incident radiation, as the resulting electron temperature is roughly equal 
to the mean photon energy of the ionizing radiation. In the calculations presented herein 
we assume a single power-law with an ene rgy spectrum L y oc v~ a between 13.6 eV and 13.6 
keV and a = 1.5 as a fiducial model (e.g. ISim et all 120081). More realistic inc ident spectral 



distributions (those used, e.g., in lEverettll2005uSimll2005tlSchurch &: Dondl2007l ) can be easily 



accommodated within XSTAR likely leading to results that are not substantially different from 
those of the next sections, although we defer such studies to future publications. 



2.2.2. Radiative Transfer in the Wind 



As noted above, equation ( 133]) is strictly valid when the attenuation of the radiation 
field by the wind itself is small; it is used in the computation of the ionization state only for 
the innermost, highly ionized regions of the wind, while at larger radii we use the absorbed 
radiation field is employing XSTAR for the computation of the relev ant opacities; the metho d 
is described below and it is very simila r to the recent work of, e.g., Schurch fc Pond ( 120071 ). 
Dorodnitsyn. Kallman &: Progal ( 120081 ) and ISchurch. Done &: Progal (120091 ) . Since at present 
we are interested mainly in the absorption properties of the above configuration, whose 
Thomson opacity is generally small, we treat the gas along the LOS independently as a 
collection of radially discrete slabs of constant (local) density, and the total column depends 
only on the LOS angle 9 in accordance with Figure |2h. At large distances where the value 
of the ionization parameter drops and He and H opacities become important, we let XSTAR 
subdivide each zone as necessary for the proper treatment of the radiative transfer; since 
in this case the absorption mean-free path at the characteristic frequencies is much smaller 
than the local radius, the plane parallel geometry employed by XSTAR for this purpose is 
adequate. 

With the ionizing luminosity L and a spectrum specified at x — 1, the entire LOS wind 
with the radial length scale of log(Ax) ~ 10 is divided up into a number of (plane parallel) 
slabs/zones (typically ~ 40) of logarithmically-equal thickness in radius; i.e. Ax/x = const. 
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Then, using XSTAR we calculate the absorption and emission coefficients within, say, the 
(i + l)-th zone, using as input the (properly normalized) output spectrum of the i-ih zone, 
Lj. We thus obtain the forwardly emitted spectrum in emission lines Lf+f 1 and in continuum 
^(cont) ^ computing the free-free, free-bound and bound-bound atomic transitions within 
the (i + l)-th zone (see XSTAR manual), while we ignore for this calculation the part of the 
backward directed flux. We also compute the fraction of the radiation flux of the i-ih zone, 
Li, that is transmitted through the (i + l)-th zone, Lf+\, given by 

Lfll = Lte-W , (34) 

where Ti + i(Nn,i) is the optical depth (including continuum and line photoabsorption as well 
as Thomson scattering) of the (i + l)-th zone, with Nji,i being the total hydrogen equivalent 
column of this slab. Finally, we produce the (properly normalized) spectrum of the (i + l)-th 
zone, by adding the forwardly emitted and transmitted spectra, i.e. 

T — 7~( tr ) _i_ r( cont ) i r( line ) (or\ 

J^i+i — J^i + i t -I- J^i + i , yoo) 

with the process being iterated with Lj + i in place of Li, to compute Lj + 2 and so on. 



2.2.3. Modeling Absorption Line Profiles 



The above procedure can be applied to a specific atomic transition to produce the re- 
sulting line profile and an absorption feature since this is the dominant process along the 
observer's LOS to the point-like ionizing source. The corresponding emission is distributed 
isotropically and makes a negligible contribution to specific intensity at the particular fre- 
quency. The radiation transfer of a specific feature is dominated by the optical depth of the 
wind at the specific frequency r(v), which is expressed by 



r(v) = a{v)N H {p) 



with the line photoabsorption cross-section 



a = 0.0U95(fij/Au D )H(a,u) 



(36) 



(37) 



where is the oscillator strength of the transition between the i—th and j—th levels of an 
ionic species and Au D is the Doppler broadening factor estimated by Av D w (At> turb /c)z/ 
relative to the centroid (rest-frame) frequency z/ with f tU rb being either thermal or turbulent 
velocity of the medium, assuming that the o yerall absorption profiles are well a pproximated 



locall y by the Voigt function H(a,u) (e.g. iMihalasI Il978l : iKotani et al.l 12000 



Hanke et al. 



20091 ) . This represents a line profile whose centroid is dominated by the Doppler broadening 
while the wings are characterized by the damping/Lorentzian profile and defined by 



H(a, u) = — 

7T 



~ y2 dy 



(u — y) 2 + a 2 



(31 
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In the above expression, we use a = T /(A.'k/S.vd) where V is the Einstein coefficient and 
u = (y — v )/Av D is the dimensionless frequency spread about the transition frequency. 
Irrespective of their thermal or turbulent velocities, the winds we consider provide a well 
defined velocity shear between two adjacent slabs, whose effect on the radiative transfer of a 
line photon is very similar to that of the turbulent velocity. To rid the situation of as much 
unwanted uncertainty as possible we then consider the thermal and turbulent broadenings 
to be negligible and we consider only the shear broadening, so that 

Av D « (Av tnTh /c)v = (Av sh /c)u , (39) 

with Af s h being the LOS velocity shear. In addition to the broadening, the velocity of the 
absorbing medium shifts also the transition frequency blueward as u' Q = u /[l — v\ os (x, 9)/c] 
where v\ os (x,6) = v r (x,9) is the flow velocity along LOS and u' is the overall observed 
frequency. With these considerations the expression for the dimensionless frequency spread 
about the specific transition reads 

U (X,9) = U/U °- 1/ } 1 - / VUX)/C] . (40) 
Av sh /c 



2.2.4- Modeling the Absorption Measure Distribution (AMD) 

In a wind model with well-defined density as a function of radius along the LOS, it 
is natural to use the value of £, the parameter that determines the presence (or not) of a 
given ion, as a proxy for the radius r along the wind in order to study the structure of its 
ionization. Then the hydrogen equivalent column of specific ions, Nh, as a function of £ 
can be used to determine the wind density distribution as a function of radius. This has 
been precisely the approach of HBK07, who defined the AMD as the hydrogen equivalent ion 
column density per logarithmic ionization parameter interval; in view of our specific wind 
models, this has the following dependences [Eqns. (1301) and (|33|) ] 

"'ri' TO ^"" 1 '^' < 41 > 

d[log£(r,0)j 

Therefore, depending on the specific run of the wind density with radius (via q), the AMD 
is constant, i.e., independent of the radius r or ionization parameter of the plasma £ when 
5 = 1, while it scales as AMD oc r -1 / 2 oc £, i.e., monotonically increasing with ionization £ 
for the BP82 models (q = 3/4). 



3. Results 



In this section we discuss in detail the implementation of the procedures described 
above and the ensuing results concerning the properties of individual atomic transitions as 
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well as their dependence on the parameters of the problem. In particular we focus on specific 
ionic charge states that are important in that they can serve as spectroscopic diagnostics 
of the ionization state of the wind medium. Because of the specific, continuous decrease 
of the MHD wind density inherent to the accretion-origin wind scenario, our models span 
a wide range of ionization parameter space and hence a large number of the corresponding 
atomic transitions (in absorption). Their global properties are encapsulated in the AMD, 
namely the distribution of (locally) absorbing column per logarithmic ionization parameter 
interval, a quantity originally introduced by HBK07, which can and has been compiled for 
a number of AGNs usin g comprehen sive fits to the ensemble of absorption features in their 



up 

entire X-ray spectrum (iBeharl I2009I ). Furthermore, given that to each ionization state of 
the plasma corresponds also a specific wind velocity, there exists a correlation between the 
presence of a specific atomic transition and the corresponding plasma velocity. Employing 
this generic property, we demonstrate a progressive absorption of the incident X-ray by 
matter outflowing at the different velocities (as discussed in §2.2), to produce the eventual 
absorption line profile of specific transitions. Finally, because of the strong dependence 
of the wind column density along the LOS angle 9, as discussed in §2.1, we consider two 
representative cases for comparison [employing always the q = 1 field geometry and the 
parameters used in producing Figs [TUB] , namely 6 = 30° (high latitude LOS) and 60° (low 
latitude LOS). 

Considering the underlying complexity of the observed broad-band SEDs of AGNs we 
have decided to adopt the simple prescription of a single power-law form for the ionizing flux 
discussed above and use different values of its slope a to simulate the diversity of observed 
AGN SEDs. This approach, while it should be viewed with caution vis-a-vis the detailed 
properties of the wind ionization, we believe that it is sufficient in capturing the trends of 
the winds' ionization properties with changes in the overall AGN SED. The photoionization 
calculations we present include the most import ant elements (H, He, Ca, Na, O, Ne, Mg, Si 



and Fe) with abundances set to the solar values (IGrevesse. Noels fc Sauvallll996l ). Based on 
this set-up the major characteristics of our hydromagnetic wind are listed in Table [1] which 
we shall discuss in subsequent sections. 



3.1. Absorption Measure Distribution (AMD) 

The unique characteristic of any wind model ionized by a point-source located at its 
origin is the radial distribution of density and therefore the corresponding distribution of 
the ionization states of the different elements seen in absorption in the source continuum. A 
particular wind model is therefore characterized by the distribution of the ionization states 
of the different elements, their velocities and the distribution of their column densities or 
their (local) equivalent hydrogen column Nh', these determine the run of the wind's density 



- 19 - 



Table 1. Characteristic Wind Properties of Our Models. 



Parameter a 


Fiducial Value (0 = 30°/60°) 


Physical Significance 


M 


10 6 


Mass of the central black hole in units of Mq 


a 


-1.5 


Power-law index of the incident spectrum (in energy) 


Lx 


3.3 X 10 42 erg s" 1 


Incident X-ray luminosity 


r a 


~ r s 


Radius of the innermost foot point of the wind 


q 


1 


Scaling of magnetic fields 


rh 


0.1 


Dimcnsionless (conserved) mass-accretion rate 


Vw 


0.5 


Outflow (wind) rate relative to mass-accretion rate m 


e 


0.2 


Radiative efficiency 


^in 


6.7/2.0 


Dimcnsionless LOS Inner radius of the wind 


7l(x in ) 


1.9 x 10 9 /4.5 x 10 10 cm- 3 


Density of the wind at x = X[ n 


Ax 


~ 10 10 


LOS length scale of wind 


logC(^in) 


8.62 / 8.3 


Ionization parameter at x = x; n 


ui os (Fexvii) 


~ 100 - 300 km s- 1 


Characteristic LOS velocity of Fe xvn 


log£(Fexvn) 


~ 2.2 - 3 


Ionization parameter of Fe XVII 


logi(Fexvn) 


~ 6 - 7/5 - 6 


LOS distance of Fe XVII 


t)los(Fexxv) 


~ 2, 000 - 4, 000/1, 000 - 3, 000 km s" 1 


Characteristic LOS velocity of Fe XXV 


log£(Fexxv) 


~ 4 - 5 


Ionization parameter of Fe XXV 


logz(Fexxv) 


~ 4- 5/3 -4 


LOS distance of Fe XXV 


"los(OVIl) 


~ 40 - 150 km s" 1 


Characteristic LOS velocity of VII 


uios(Ovm) 


~ 100 - 600/100 - 300 km s" 1 


Characteristic LOS velocity of VIII 


AN H 


~ 2.6 x 10 21 /1.8 x 10 22 cm- 2 


Local wind column density per slab 


N h 


~ 3.9 x 10 22 /2.5 x 10 23 cm" 2 


Integrated wind column density over — 1 < log £ < 4 



a See the text for the template parameter values in detail. 
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Fig. 4. — Simulated hydrogen-equivalent (local) column density (per slab) AiV#(ion) for 
major ionic elements (N, O, Ne, Mg, Si and Fe) and wind temperature T (right axis) in the 
fiducial model for (a) 30° and (b) 60°. Analytic self-similar (SS) wind temperature is also 
shown. Corresponding LOS (normalized) distance x(= r/rs) is also shown in the upper axis. 
The parameters are the same as in Figure [TJ Rapid decrease of ANh in the highest (lowest) 
end of £ is due to fully stripped ions (neutral ions). [See the electronic edition of the Journal 
for a color version of this figure.} 



and velocity with the radius r. For example, in the accelerating phase of a spherical wind, 
both the ionization parameter and the velocity increase with distance and "freeze" once it 
has achieved its terminal velocity, while the column density of ions decreases faster than 1/r 
if an ion is present in the wind accelerating stage, leveling off to iV# oc 1/r once terminal 
velocity is achieved. 

The characteristic properties of the winds presented in §2 are drastically different from 
those of spherical winds. In particular: (i) The ionization parameter £ and the corresponding 
LOS velocity t>i os decrease with distance as r -1 and r _1//2 respectively to encompass a large 
number of ionization states of many elements at different velocities, (ii) The specific column 
density per decade of radius translates for any element to a (local) equivalent hydrogen 
column ANh that is independent of the radius r. 

For a given wind model and ionizing spectrum, one can produce the corresponding 
AMD profiles, i.e. the global distribution of its ionization properties which can then be 
compared to observations. We do so in this section and compare the results to the recent 
AMD analysis of IRAS 13349+2438 observations (HBK07). The results of our calculations 
are shown in Figure H] where we present the AMD for the q = 1 wind models, with the 
ANh plotted as a function of log£ for 9 = 30° in (a) and 60° in (b). Also superimposed 
is the plasma (gas) temperature logT computed by XSTAR under local thermal equilibrium 
conditions along with the analytic self-similar (SS) temperature (oc 1/r). The ionized wind 



-21 - 



temperature decreases monotonically from T ~ 10 6 K down to T ~ 10 3 ' 5 K as £ decreases 
from log£ ~ 5 to log£ ~ -1 a range relevant to the observable ionic states of the major 
elements. For reference, the corresponding (dimensionless) LOS distance x = r/r$ from the 
central engine is also shown in the upper x-axis where in this run rg = 3 x 10 11 cm. As seen, 
the local column densities of the irradiated ions (shown here are the six major elements: N, 
O, Ne, Mg, Si and Fe) are found to be all distributed at constant value of ANh over many 
decades of £ (and therefore a similar range in LOS distance too), namely ANh ~ 2.6 x 10 21 
cm -2 for 30° and ~ 1.8 x 10 22 cm -2 for 60°, reflecting the different wind density at these 
different LOS (see Fig. |2^i). This results in an integrated column of Nh ~ 3.9 x 10 22 cm -2 
for 30° and ~ 2.5 x 10 23 cm -2 for 60° over — 1 < log£ < 4. The constant value of AMD is a 
unique characteristic of the q = 1 models with the value corresponding to case (a) (6 = 30°) 
in good agreement with the observed AMD of IRAS 13349+2438 analyzed by HBK07. 

As shown here, the computed temperature has a radial dependence (gradient) close to 
the desired 1/x one for x ^ 10 5 , but it is constant for x & 10 5 . However, because the 
thermal pressure is smaller than the wind ram pressure over the entire domain, the thermal 
pressure effects on the wind density profiles are not significant. We have also compared the 
heating/ cooling time scales of XSTAR and were found to be shorter than the local dynamical 
time scale r/v, indicating that the temperature distributions of Figure H] are indeed correct, 
even though not dynamically important. We have obtained by another calculation a wind 
solution similar to Figure |U in the cold- flow limit (K = 0). 

It should be noted that the decrease of ANh at high values of £ is due to the fact that 
XSTAR by default does not provide the column density of fully ionized (bare) ions. This 
apparent decrease is just due to the fact that all considered elements become progressively 
(completely) ionized at high £ values. At the other end of the plot, i.e. for small values 
of £, AN H appears to decrease too because the wind becomes too cold (gas temperature 
of T < 10 3 ' 5 K) to participate in radiative processes, i.e. each species turns neutral, an 
ionization state also not included in the calculation of the hydrogen equivalent column. The 
calculated AMD is otherwise constant, as expected. The fact that the corresponding ANh is 
the same in this regime irrespective of the elements used in the computation is a testimony 
that our models implement XSTAR correctly, so that the sum of the abundances of all ions 
of any element, when corrected for the elemental abundance, yields the same value of ANh- 
We also see that the (relatively) lighter elements (say, oxygen) become fully ionized at lower 
£ than heavier ones (say, iron), also as expected. In essence, the overall trend of the model 
AMD in both LOS cases (a) and (b) appears to be very similar to one another (i.e. constant 
value of ANh independent of ionization state) except for the fact that the local column ANh 
for 30° in (a) is about a factor of ~ 10 lower than that for 60° in (b) because of the angular 
dependence of the wind density function M{9) (also see Fig. [2^). 

It is also of interest to show the differential form of Figure HJ i.e. the columns of specific 
charge states of a given ionic element AN lon . This is shown in Figure |5] for iron again 
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Fig. 5. — Simulated local iron column density ANp e distribution for various iron charge levels 
(Fein, Fe VI, Fex, Fexill, Fexvn, Fexx, Fexxv, and Fexxvi) and LOS velocities v\ os (right 
axis) for (a) 30° and (b) 60° with the same fiducial model parameters as in Figure HI [See 
the electronic edition of the Journal for a color version of this figure.} 



for two different LOS. This figure shows the emergence of various ionization states of iron 
(shown here are Fe ill, Fe VI, Fe x, Fe xm, Fe xvn, Fe XX, Fe xxv and Fe xxvi) along with the 
corresponding LOS velocity v\ os for 9 = 30° in (a) and 60° in (b) corresponding to Figure HI 
As discussed above, the higher ionization states are present at the inner sections of the 
wind (smaller x), being replaced by lower ionization ions at larger radii with their columns 
given in absolute values ANp e instead of the hydrogen-equivalent values ANh- One should 
note that the maximum value (peak) of ANp c for a specific charge state is approximately 
cons tant until one gets to th e lowest ionization states, located at the largest distances (see, 
e.g., Murray fe Chiang! Il998l . for a similar transition) because the fraction of neutral iron 
relative to the total iron starts increasing as mentioned earlier for Figure HI Of interest is 
also the LOS wind velocity v\ os that corresponds to a given value of £. While t>i os decreases 
like oc £ 1//2 at high values of £ (as £ oc r _1 and v\ os oc r" 1 / 2 ), it deviates from this dependence 
at the largest radii because of the additional reduction of the X-ray flux due to absorption. 
The velocity at which each ion has its maximum column and the range in v\ os over which it 
has substantial column are of interest because they determine the absorption line profile of 
the specific ion. Thus, in the case of (a) 9 = 30°, higher charge states such as Fexxv are 
found to be blueshifted by a relatively high outflow velocity, v\ os ~ 2, 000 — 4, 000 km s _1 
at log £ ~ 4 — 5 (log x ~ 3 — 4) while lower charge states such as Fe xvn are blueshifted 
by only t> los ~ 100 — 300 km s -1 at log£ ~ 2.2 — 3 (logx ~ 6 — 7). We emphasize that the 
latter (Fexvil) line velocity and column density are also in good agreement with those of 
IRAS 13349+2438 (HBK07). For 9 = 60° in comparison, the column density increases ~ 10 
times but the corresponding line velocity and ionization stage remains roughly the same (but 
at smaller distances) because both the poloidal velocity and its projection to the LOS are 
smaller at lower latitudes (see Fig. [2b). 
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Fig. 6. — Same as Figure |5] but for oxygen. [See the electronic edition of the Journal for a 
color version of this figure] 



To conclude this section, we show in Figure [6] the column densities of specific charge 
states of oxygen again for two different LOS at (a) 9 = 30° and (b) 9 = 60°, along with the 
corresponding velocities overplotted as in Figure [5j As expected the corresponding outflow 
velocities are smaller for O VII (v ~ 40 — 150 km s _1 at log£ ~ 1 — 2 or logx ~ 7 — 8) and 
larger O Vlll (v ~ 100 - 600 km s" 1 at log£ ~ 6 - 7 or logo; ~ 2 - 3) for (a) 9 = 30°, while 
similarly for (b) 60° we find v ~ 40 — 150 km s _1 at log£ ~ 2 (log a; ~ 6) for Ovil and 
v ~ 100 — 300 km s _1 at log£ ~ 2.5 — 3 (logx ~ 5 — 5.5) for O Vlll, with the former having 
lower column per ion than the latter, as expected. 

One should note that the precise correlations between the wind variables and the ob- 
servables (i.e. AiV ion , t>i os ) in the present model are primarily a function of the specified 
mass-accretion rate m and LOS angle 9 as well as the incident X-ray spectral index a and 
therefore the model has the freedom to alter the resulting AMD distributions depending on 
the intrinsic properties of AGNs (here we focused on IRAS 13349+2438 clS db CclSG study). 



3.2. Absorption Line Spectra 

As discussed in the previous subsection and exhibited plainly in Figures [5] and EJ spe- 
cific charge states of the different elements are present at different values of £ to which 
correspond different values of column density and velocity. We have also indicated that the 
simplest assumption about the optical depth of a specific transition is that it is determined 
by the velocity shear of the wind, a quantity that is also calculable within our models [see 
equation (37)]. Armed with this information we present in this section the profiles of spe- 
cific atomic transitions in absorption. With the spectral emissivity and velocity given, we 
are also in a position to calculate the corresponding emission profiles. However, since the 
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Fig. 7. — Simulated broad-band X-ray spectra L v for (a) 30° and (b) 60° as a sequence 
of progressive X-ray illumination: an initial power-law (slab 1) and the subsequent spectra 
(as indicated by ionization parameter log£-value of each slab) with the same fiducial model 
parameter as in Figure H] [See the electronic edition of the Journal for a color version of this 
figure.] 

most prominent AGN lines lie in the O-UV part of the spectrum, a comprehensive approach 
would require also a more precise model of the UV part of the spectrum; as such we defer 
line emission to a future work. 

A demonstration of the progressive effect of absorption on the transmitted spectra is 
presented in Figure [7J This figure exhibits the spectra between 1 eV and 10 keV of the 
radiation transmitted past the innermost slab 1 (power-law) through the subsequent slabs 
indicated by the values of log£ for two different LOS angles, (a) 9 = 30° and (b) 60°. Ab- 
sorption by the consistently ionized wind becomes apparent first at energies E <^ 1 keV, 
becoming progressively deeper and shifting to lower energies as the radiation propagates 
through larger wind column. Absorption also increases with the inclination angle 9 as one 
would qualitatively expect. Note that in this simplistic model the incident radiation includes 
only the X-ray power-law continuum ignoring all the other (potentially rather important) 
components. Also, it ignores emission that eventually may come into our LOS due to reflec- 
tion and scattering in the wind. For this reason, the deep absorption imprinted in Figure [7J 
should not be interpreted literally in a quantitative sense, but as a qualitative indicator of 
the absorption function of £ along a given LOS. 

Let us turn our attention to specific atomic transitions, namely those of Fexvil and 
Fe xxv of column AiVpe with velocity function of £ as given in Figure |5] and shear 

given by our model. Using the expressions for the photoionization cross-section and Voigt 
function [eqns. (|36|) - (158]) ] one can construct the characteristic absorption profiles; these 
include both the Doppler blueshift due to LOS velocity v\ os and the line broadening of the 
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Fig. 8. — A sequence of the simulated absorption features (indicated by ionization parameter 
log c;-value of each slab) for Fe xvn (upper panels) and Fe xxv (lower panels) as the incident 
X-ray (power-law shown by dotted horizontal lines) transmits through the LOS wind of 30° 
(left panels) and 60° (right panels) with the same fiducial model parameters as in Figure HI 
Vertical dashed lines denote the rest-frame wavelength. 

Voigt function in equation fl38|) . Figure [8] shows a sequence of the simulated absorption line 
for Fexvil with the rest-frame wavelength at 15.014A (upper panels) and Fexxv with the 
rest-frame wavelength at 1.850A (lower panels) for two LOS angles, 9 = 30° (left panels) and 
60° (right panels). We focus here on these prominent iron charge states because they are less 
contaminated by other adjacent spectral features. We adopted the following values for the 
relevant atomic data: / y = 0.69/2.3 and V = 4.48 x 10 14 /2.27 x 10 13 s" 1 for Fexxv/Fexvil, 
respectively. The figure provides the spectra transmitted through increasing wind column 
along the LOS. The absorption line obtains its largest blueshift and broadest contribution by 
ions at the higher values of £, becoming progressively deeper, less blue-shifted and narrower 
as the contribution to the absorption is affected by regions of lower £ (larger r) which provide 
the highest column of the specific ion (and also smaller velocities and shear). Again, the 
spectra are broader and deeper for the higher LOS angles reflecting higher (local) velocities 
at specific values of £ and higher overall columns. 

In concluding this section we would like to point out the work of Gabel et al. (2003) 



-26 - 



on the UV absorption features of NGC 3783 which were found to have velocity structure in 
agreement with that of the corresponding X-ray ones, thereby arguing for the consistency of 
the entire wind structure. Also, Collinge et al. (2001) have found that the low ionization Fe 
absorption features of NGC 4051 had corresponding UV absorption features, while the high 
ionization, higher velocity X-ray absorption features of the spectrum lacked an equivalent 
UV absorption, indicating the absence of these ions in the higher ionization, higher velocity 
plasma. 



4. Summary and Discussion 

In this work we have presented a detailed study of the ionization structure of model MHD 
winds off accretion disks; this is our first attempt to model within this context the recent 
observations of absorption features in the ASCA, XMM and Chandra X-ray spectra, the so- 
called warm absorbers. To this end we have employed the self-similar 2D hydromagnetic wind 
models developed by CL94 that provide the fluid 2D density and velocity fields, which we 
coupled to photoionization calculations using XSTAR. Our consideration of magnetocentrifugal 
(MHD) winds/outflows in this work has been motivated and supported in part by recent 
observational implications that at least in a number of AGNs and GBHCs the inferred driving 
mechanism of the observed X-ray ionized wind medium is magnetic rather than thermal or 
radiative (e.g., see Miller et al. 2006, 2008 for GRO J1655-40; Kraemer et al. 2005 and 
Crenshaw & Kraemer 2007 for NGC 4151). 

Given the scope of our paper, our model wind is necessarily overly simple. It ignores a 
host of issues that affect the structure of winds off accretion disks and replaces them with 
the self-similar models of CL94. The interested reader can get a feeling of the multitude and 
com plexity of the issue s not addressed in t he pr esent treatment by looking at the works of, 



S3 » and bwa fc Ka^anl (EJ (and reference, therein) and also 



(120091 . for disks with several distinct values of the accretion rate m). These include, among 



others, the radiative transfer in the wind, the ensuing effects of radiation pressure and also 
the influence of the vertical gradient of B\ in the launching of the wind. As shown in these 
references the winds can be radiation driven in the inner part of the disk (due to the increased 
flux at this spatial domain) while being magnetically launched at larger radii. In this respect 
o ur model wind s are quite different in structure and variability from the more realistic winds 



of lProgal (120031 ) at small radii, while they should be more similar to those of his at larger radii 
where radiation pressure is smaller and the wind is magnetically launched. Realistic models 
must by necessity consider the transition of the underlying disk from radiation pressure, 
at small r, to gas pressure dominance at larger radii; this fact would likely force a different 
choice of our boundary condition -i//(90°) and therefore yield a field geometry that breaks the 
self-similarity of our solutions. Wind models very similar to those used herein but with the 
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inclusion of t he effects of the radiation pressure that break the self-similarity were discussed 



by KK94 and lEverettl (1200a ) 



One should therefore view our models with the above caveats in mind. The significance 
and assumptions of the models we propose is basically justified a posteriori by their ability 
to interpret the observations. Self-similarity is one of the fundamental assumptions of our 
model winds. However, whether exactly self-similar or not , any model t hat would attempt 



to account the entire range of ionic species shown e.g. in iBeharl ( 120091 ) must by necessity 
cover a very broad range in the photoionization parameter £. Then, the column of each such 
ion provides a measure of the corresponding hydrogen equivalent column iV# as a function 
of £, i.e. the AMD. These two quantities (Ng,£) can then be employed to provide a measure 
of the gas density n(r) as a function of the distance r from the X-ray source. Barring 
the possibility of several independent regions at different distances but similar columns as 
indicated by the functional form of the AMD (one could consider this is possibility at the risk 
of introducing an inordinate number of free parameters), models with radial density profiles 
as those considered here are a natural consequence of the AMD obtained by the Chandra 
observations. 

The wind ionization structure was followed in ID as described in §2.2, namely along the 
observer's LOS, assuming the ionizing source to be point-like. Clearly, a more comprehensive 
treatment of this problem should take into account also the backward emitted radiation which 
will reach the observer from the regions of the disk on the other side of the black hole, as 
well as the scattered radiation, whose effects could be significant. We plan to return to these 
issues in a future publication. 

The 2D structure of the winds considered in our analysis implies, at a minimum, a 
two-parameter description of the absorption features. However, the self-similarity of the 
problem simplifies further the treatment by allowing the separation of the r and 9 variables 
with the wind density obtaining the form n(r,9) oc r 2q ~~ 3 Af(8). The parameter q is a most 
important parameter of these models because it determines both the radial dependence of 
the ionization parameter £ and the wind column density per decade of radius along a given 
direction (LOS). In the models examined in the present work we restricted our study to 
the value q — 1, which leads to a radial density profile n(r) oc 1/r. This is an interesting 
profile in that it provides for an ionization parameter £ with a similar dependence, i.e. 
£ oc 1/r, and most importantly, with equal column per decade in radius, across the entire 
range in radius and the corresponding range in £; as a result, ionic species of very different 
ionization properties, existing over a wide range of ionization parameter (and radius), have 
roughly comparable column densities, independent of the distance and are therefore possible 
to detect. As noted in CL94, the value of q determines also the axial current distribution 
in radius within the wind, with the corresponding magnetic energy per unit length at the 
wind base being also constant per (cylindrical) decade in radius. It is worth comparing our 
work with that of KK94, who used the same type of MHD wind with the value of q as we do 
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herein: while these authors focused their study on the effects of dust on radiation transfer, 
we focus our attention on the effects of the ionizing radiation on the X-ray spectra of AGNs, 
in particular on the ionized absorbers. KK94 also noted in passing the relevance of their 
models to the X-ray AGN spectra well ahead of the detailed outflow observations made with 
Chandra. 

While the structure of the wind in r and 9 is provided by the models of CL94, for 
the remaining parameters of our models, namely the normalization of the wind density and 
the X-ray luminosity of the ionizing source we have chosen to use the scalings of ADAFs, 
or better ADIOS; as a result the X-ray luminosity is proportional to the square of the 
(normalized) accretion rate at the inner edge of the disk m 2 rather than simply m, assumed 
in standard accretion disks. Therefore, restricting ourselves to models with q = 1, the global 
ionization structure of a given wind, including the normalization of the column density, 
depends only on two parameters, namely m and the observer's inclination angle 9. This 
parametrization provides an extremely economical set of assumptions concerning not only 
the outflows, i.e. the wind ionization structure seen in absorption in detailed X-ray spectra, 
but also for the entire (radio-quiet) AGN unification picture: Figure [2] which exhibits the 
normalization of the wind density as a function of 9 along the magnetic field line of \l/ = \l/o, 
makes apparent the difference in column between face-on and edge-on views, implying that 
the wind, if extending to sufficiently large radii, can in fact serve as the proposed molecular 
torus associated with AGN unification, a point originally made by KK94, also for winds with 
q = 1; it is of interest to note that the few ob jects for whic h sufficiently detailed observations 
exit are consistent values q ~ 1 (HBK07; iBeharl 120091 ). Within this same context and 



ionizing luminosity considerations, one should note the dependence of source obscuration 
on the X-ray luminosity, namely its reduced value for objects accreting at a higher fraction 
of their Eddington rate rh, which apparently is in general agreement with observations. It 
remains to be seen whether these notions can withstand the scrutiny of more consistent and 
encompassing observational tests. 

At this point we would like to stress the importance of the AMD in the study of AGN 
outflows/winds, a quantity enunciated by HBK07 and modeled in detail in the preceding 



sections. Analyses similar to those of HBK07 and IBeharl (120091 ) are indispensable because 
they produce a consistent analysis of the entire set of absorption features in the AGN X-ray 
absorption spectra. At the same time they underscore the unique value of X-ray spectroscopy 
which, in a wavelength band of ~ 1.5 decades, can encapsulate the properties of ions that 
span ~ 5 orders of magnitude in ionization parameter and, for models with q ~ 1, a similar 
range in radius; at the same time, measurement of their absorption columns yields the 
equivalent hydrogen column of the flow, Nh, over a similar range in radius, thereby going 
a long way toward the determination of the physics underlying the outflow dynamics. The 
fact that the AMD analyses to date are roughly independent of £ provide support to our 
use of the q = 1 or n(r) oc 1/r models, reiterating that a similar density profile was invoked 
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on the basis of AGN and GBH C timing properties (IPapadakis. Nandra &: Kazanad 12001 
Kazanas. Hua fc Titarchukill996l ). 



We have thus presented a concrete example of the AMD dependence on £ for our q = 1 
models with e = 0.2, rh = 0.1 and two different values of the observer's inclination angle (see 
tabled] for details). As expected we found the AMD to be constant, i.e. independent of £, 
over many decades in this parameter (i.e. the local column density ANh is independent of 
ionization states £ of ions). With rh = 0.1, rjw — 0.5, equation f )30|) implies AN H ~ 2.6 x 10 21 
cm" 2 (yielding a total column of N H ~ 3.9 x 10 22 cm" 2 ) for 9 = 30° and AN H ~ 1.8 x 10 22 
cm -2 (total of Nh ~ 2.5 x 10 23 cm -2 ) for 60° over —1 < log£ < 4 erg cm s~ \ in good 
agreement with the observed AMD of IRAS 13349+2438 (HBK07; iBeharl l20Q9h and also 
some other AGNs detailed below. One should note that these values depend primarily on 
rh and 9 and are independent of the black hole mass. The mass of the object gets involved 
only as a measure of its total luminosity, which does not appear in the expression for Nh, 
implying that these models could in principle be applicable also in GBHCs. High quality 
X-ray absorption data are in fact available for the GBHC GRO J1655-40 and were used to 
argue for magnetic driving of the wind in this system too (IMiller et al.ll2006l . 120081 1 . These 
spectra are distinguished from those of AGNs by the prominent absence of low ionization 
state ions. This is to be expected given that the presence of the companion star limits the 



12 



cm. 



extent of the disk to roughly half the distance between the two objects or r ~ 10 
Given that the Schwarzschild radius of the compact object is r$ — 10 6 cm the entire disk 
size covers only a range of log x ~ 6 in radius; considering (based on Fig. H]) that in the inner 
region of the wind (r ~ 1000 — 3000r,s ~ 10 9 cm) major elements (even heavier species) are 
almost fully ionized, one would expect the presence of ions over only a factor of 1000 in £, 
in rough agreement with the observation. 

Given that our models provide also the complete velocity field of the MHD wind we have 
also produced a sequence of synthetic absorption profiles as described in §2.2 and shown in 
§3.2. We have done so for two charge levels of Fe, namely Fexxv and Fexvil. With our 
fiducial model predicting the Fexvil column to be maximum at log£ ~ 2.2 — 3 erg cm s _1 , 
we then infer the corresponding LOS velocity (see Fig. [5]) to be t>i os ~ 100 — 300 km s" 1 in 
excellent agreement with the observed values of IRAS 13349+2438. It is important to note 
that the model produces not only the correct velocity v\ os at the maximum value of AAp e (£) 
but also the observed ANh{C) normalization for the same inclination angle. It is therefore 
possible with measurements of the combined line widths and column densities to reproduce 
within the present models both q, 9 and rh, thereby providing a complete specification of 
these winds. In addition to IRAS 13349+2438, the nearby bright Seyfert MCG 6-30-15 has 
been observed to show, as one of the multiple ionization zones, an X-ray ionized absorber 
with an outflow v e locity of ~ 1900 km s~^ at 1 ok£ ~ 3.85 and N H ~ 9 x 10 22 cm -2 (e.g. 
Young et all 1200.4 iHolczer. Behar &: Aravl 120091 ) also in good agreement with our model 
results for rh ~ 0.1 (see Fig. |5]). Clearly, a different choice of a set of conserved quantities 
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(wind variables) in the model would produce a slightly different field line geometry in which 
the resulting MHD outflow could in principle obtain higher LOS velocities (perhaps by factors 
of magnitudes), and this will be studied in detail in a future work. 

As noted above our model is quite success ful in reproducing the observed AMD of 
IRAS 13349+2438. More recently, iBeharl d2009h presented a compilation of the AMD of 
a number of AGN with dNn/dlog^ which are slightly different from constant but with 
q— values of n(r) oc r 2q ~ 3 which are still very close to unity; e.g. NGC 3783 (q ~ 0.89), 
NGC 5548 (q ~ 0.94), NGC 7469 (q ~ 0.9) and MCG-6-30-15 (q ~ 0.95). We believe 
some of these AMDs are sufficiently close to those of our fiducial model to be virtually 
indistinguishable. By comparison, the MHD wind models of BP82 have q = 0.75 and 
n{r) oc r~ 3 / 2 . Their ionization parameter decreases more slowly with r (£ oc r -1 / 2 ) and hence 
a range of 5-6 decades in £ implies a range of 10-12 decades in radius with the corresponding 
distances being unrealistic. Most importantly, the wind columns would decrease with radius 
(Nh oc r~ 1//2 ) and the AMD dependence on £ would be dNjj/dlog^ oc £, in clear disa greement 



with these observations, and thus their model is essentially ruled out (iBeharl 120091 ) 



Despite the apparent success of this first model, one should bear in mind that several 
aspects have been treated in a rather simplified fashion. Here we discuss some of them and 
their influence on the results presented so far: 

1. The shape of the ionizing spectrum used so far (F v oc v~ a between 13.6 eV and 
13.6 keV with a = 1.5) i s quite simplisti c, however, such spectra are often used in simi- 



lar type calculations (e.g. ISim et al.ll2008l ). This is a crude approximation to the observed 



spectra characterized by complex s pectral shap es. Spectral features like the BBB and the 
UV to X-ray luminosity ratio «ox ( jElvia 120001 ) play a role in our results (a recent review 



of the broad band AGN SEDs can be found, e.g., in Risaliti Sz Elvis 20041 ). More recently. 



Grupe et al.l (120041 ) sampled 110 bright soft-Xray selected AGNs for a simultaneous study 
of Optical/UV and X-ray data and reported strong correlations between the X-ray spectral 
slope and the Optical/UV slope which should be included in a more compreh ensive treat- 
men t. Multi-components o f the broad-band AGN spectra, as those discussed by lElvisI (120001 ) 
and lRisaliti &: Elvisl ( 120041 ). should play an important role in characterizing the observed X- 
ray spectral features. Adding softer photons (responsible for the Optical/UV components) 
could impact the radiative t ransfer between those m ore complex photon distributions and 
the ionized matter (see, e.g.. lEverettll2005l ; ISimll2005l for mult i- component injected spectra). 



In fact, we have mimicked the presence of such softer photons by choosing larger values of the 
spectral index of our models (a > 2). For the same total luminosity (implemented in the def- 
inition of £), it was found in this trial run that the radial position of the peak column density 
of a given ion (e.g. Fexvil) decreases, because the presence of a given ion requires a certain 
flux of ionizing photons per atom; in a steep spectrum the proper ratio is found at larger 
values of the parameter £ or correspondingly smaller values of the distance x = r/r$- As a 
result, the widths of the corresponding transitions should be larger and therefore the proper 
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ionizing spectrum is necessary for the corr ect interpretation of the relations between ionic 
column densities and wind kinematics (see 



Everett 


2005; 


Sim 


2005 



2. The ionizing source has been considered, for simplicity, to be point-like. A source 
of finite size will provide for more complex illumination, especially for parts of the wind at 
which the source angular extent is significant. Most importantly an extended source would 
impact the absorption line profiles, since a line of sight may pass through matter other than 
that used in section 3.2. Finally, to treat the radiative transfer correctly at all energies, one 
should eventually need to resort to a 2D approach rather than the ID model (radial only) 
presently used, which would modify the global ionization structure (Fig. [3b). 

3. The wind equations, as presently implemented include only MHD forces, while it 
is apparent that in the presence of the ionic species we produce one should also include 
the effects of radiation pressure. Because the latter depends on the ionization state of the 
plasma which in turn depends on its kinematic s, implemented correctly, this should be done 
iteratively in a way similar to lEverettl (120051 ) . Based on the rather small effect that the 
radiation force seems to have on the structure of these winds, we believe that our present 
results are generally valid. 

In summary, we have presented above a first attempt at interfacing theoretical models 
of MHD winds with AGN observations, in particular the absorption features in their X- 
ray spectra of ionized outflows. We have placed most of our emphasis on modeling the 
corresponding AMDs, which, based on the existing observations seems to favor a specific 
value of the model parameter q which determines the radial density dependence of our models, 
namely q ~ 1. With the value of this parameter set by observation, our models present a two- 
parameter family of AGN wind structure, namely rh and 9. We have found that, based on the 
limited number of objects discussed above, our models fare rather well in accommodating the 
observations and the possibility of incorporating AGN unification within the same models, 
as discussed in KK94. There is still a multitude of issues related to X-ray absorption that 
remain open, e.g. whether the model can accommodate the h igh velocity (v/c ~ . 1 — .25) 



outflows associated with X-ray absorption features (see, e.g., ISiml 120051 ; ISim et al.ll2008l for 
some attempts with Monte Carlo simulations) in the spectra of some bright quasars such as 
APM 08279+5255, PG 1 211+143 and PDS 456 (e.g. Idiartas et al.ll2009t bounds fc Re"e~^ 



2009 



Reeves et al.l 120091 ) . In the context of our model this point may be well explained 
with an optimized magnetic field geometry [i.e. \l/(r, 9) distribution determined by GS- 
equation (|27|) ] that provides a favorable LOS velocity component. Important as they are, 
these go beyond the scope of the present work and we expect to return to them in future 
publications. 



We are grateful to our anonymous referee for his/her constructive comments that im- 
proved the manuscript. We would like to thank Tim Kallman for his help with XSTAR 



-32 - 



incisive comments. We express our gratitude to George Chartas for his comments on the 
model, Takanori Sakamoto and Javier Garcia for their assistance with scripting and running 
XSTAR as well as helpful comments. This work was supported in part by NASA ADP grant. 

REFERENCES 

Antonucci, R. R. J. & Miller, J. S. 1995, ApJ, 297, 621 

Behar, E. et al. 2003, ApJ, 598, 232 

Behar, E., 2009, ApJ, 703, 1346 

Boroson, T. A. 2002, ApJ, 565, 78 

Boroson, T. A. & Green, R. F. 1992, ApJS, 80, 109 

Blandford, R. D & Payne, D. G. 1982, MNRAS, 199, 883 (BP82) 

Blandford, R. D & Begelman, M. C. 1999, MNRAS, 303, LI 

Brandt, W. N., Chartas, G., Gallagher, S. C, Gibson, R. R., & Miller, B. P. 2009, 
la7Xiv:0909.0958V l 

Chartas, G., Saez, C, Brandt, W. N., Giustini, M., & Garmire, G. P. 2009, ApJ, 706, 644 

Contopoulos, J., & Lovelace, R. V. E. 1994, ApJ, 429, 139 (CL94) 

Crenshaw, D. M., & Kraemer S. B. 2007, ApJ, 659, 250 

Dorodnitsyn, A., Kallman, T., & Proga, D. 2008, ApJ, 675, L5 

Ferrarese, L. & Merrit, D. 2000, ApJ, 539, L5 

Elvis, M. 2000, ApJ, 545, 63 

Everett, J. E. 2005, ApJ, 631, 689 

Grevesse, N., Noels, A., & Sauval, A. 1996, in "Cosmic Abundances" ASP Conference Series, 
99, S. Holt and G. Sonneborn, eds. 

Grupe, D. et al. 2004, AJ, 127, 1799 

Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 

Hanke, M. et al. 2009, ApJ, 690, 330 



-33 - 



Holczer, T., Behar, E. & Kaspi, S. 2007, ApJ, 663, 799 (HBK07) 
Holczer, T., Behar, E. & Arav, N. 2009, (submitted to ApJ) 
Hopkins, P. F. et al. 2005, ApJ, 630, 705 

Hua, X.-M., Kazanas, D. & Titarchuk, L. G. 1997, ApJ, 482, L57 
Kazanas, D., Hua, X.-M. & Titarchuk, L. G. 1996, ApJ, 480, 735 
Kallman, T., k Bautista, M. 2001, ApJS, 133, 221 
Kawanaka, N., Kato, Y., & Mineshige, S. 2008, PASJ, 60, 399 
Kongl, A. & Kartje, J. F. 1994, ApJ, 434, 446 (KK94) 
Kotani, T. et al. 2000, ApJ, 539, 413 
Kraemer, S. B. et al. 2005, ApJ, 633, 693 
Krolik, J. H. & Begelman, M. C. 1988, ApJ, 329, 702 
Lynden-Bell, D. 1969, Nature, 223, 690 

Mihalas, D. 1978, "Stellar Atmospheres", W. H. Freeman & Company, New York, NY, 2nd 
edition, p. 279 

Miller, J. M. et al. 2006, Nature, 441, 953 

Miller, J. M., Raymond, J., Reynolds, C. S., Fabian, A. C, Kallman, T. R., & Homan, J. 
2008, ApJ, 680, 1359 

Murray, N., & Chiang, J. 1998, ApJ, 494, 125 

Nandra, K. et al. 1998, ApJ, 505, 594 

Narayan, R., & McClintock, J. E. 2008, New Astron. Rev., 51, 733 
Narayan, R. & Yi, 1994, ApJ, 428, L13 
Narayan, R. & Yi, 1995a, ApJ, 444, 231 
Narayan, R. & Yi, 1995b, ApJ, 454, 710 

Ohsuga, K., Mineshige, S., Mori, M., & Kato, Y. 2009, PASJ, 61, L61 
Papadakis, I. E., Nanadra, K. & Kazanas, D. 2001, ApJ, 544, L133 



-34 - 



Pounds, K. A. & Reeves, J. N. 2009, MNRAS, 397, 249 

Proga, D. 2003, ApJ, 585, 406 

Proga, D. & Kallman, T. R. 2004, ApJ, 616, 688 

Reeves, J. N., O'Brien, P. T., Braito, V., Behar, E., Miller, L., Turner, T. J., Fabian, A. C, 
Kaspi, S., Mushotzky, R., & Ward, M. 2009, ApJ, 701, 493 

Reynolds, C. S. & Fabian, A. C. 1995, MNRAS, 273, 1167 

Risaliti, G. & Elvis, M. 2004, in "Supermassive Black Holes in the Distant Universe", Ed. 
A. J. Barger, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2004, p. 187 

Sanders, D. B. et al. 1989, ApJ, 347, 29 

Schurch, N. J. & Done, C. 2007, MNRAS, 381, 1413 

Schurch, N. J., Done, C. & Proga, D. 2009, ApJ, 694, 1 

Sim, S. A. 2005, MNRAS, 356, 531 

Sim, S. A., Long, K. S., Miller, L., & Turner, T. J. 2008, MNRAS, 388, 611 
Tueller, J. et al. 2008, ApJ, 681, 113 

Ueda, Y., Akiyama, M., Ohta, K., & Miyaji, T. 2003, ApJ, 598, 886 
Urry, M. C. & Padovani, P. 1995, PASP, 107, 803 

Young, A. J., Lee, J. C, Fabian, A. C, Reynolds, C. S., Gibson, R. R., & Canizares, C. R. 
2005, ApJ, 631, 733 



This preprint was prepared with the A AS IAT^X macros v5.2. 



